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Abstract —This paper develops a comprehensive analytical 
framework with foundations in stochastic geometry to char¬ 
acterize the performance of cluster-centric content placement 
in a cache-enahled device-to-device (D2D) network. Different 
from device-centric content placement, cluster-centric placement 
focuses on placing content in each cluster such that the collective 
performance of all the devices in each cluster is optimized. 
Modeling the locations of the devices by a Poisson cluster process, 
we define and analyze the performance for three general cases: (1) 
fc-Tx case: receiver of interest is chosen uniformly at random in 
a cluster and its content of interest is available at the closest 
device to the cluster center, (ii) €-Rx case: receiver of interest 
is the closest device to the cluster center and its content of 
interest is available at a device chosen uniformly at random from 
the same cluster, and (lii) baseline case: the receiver of interest is 
chosen uniformly at random in a cluster and its content of interest 
is available at a device chosen independently and uniformly at 
random from the same cluster. Easy-to-use expressions for the 
key performance metrics, such as coverage probability and area 
spectral efficiency (ASE) of the whole network, are derived for 
all three cases. Our analysis concretely demonstrates significant 
improvement in the network performance when the device on 
which content is cached or device requesting content from cache 
is biased to lie closer to the cluster center compared to baseline 
case. Based on this insight, we develop and analyze a new 
generative model for cluster-centric D2D networks that allows 
to study the effect of intra-cluster Interfering devices that are 
more likely to He closer to the cluster center. 

Index Terms —D2D caching, cluster-centric content placement, 
clustered D2D network, Thomas cluster process, stochastic ge¬ 
ometry. 

1. Introduction 

D riven by the increasing mobile data traffic, cellular 
networks are undergoing unprecedented paradigm shift 
in the way data is delivered to the mobile users [2]. A key 
component of this shift is device-to-device (D2D) commu¬ 
nication in which proximate devices can deliver content on 
demand to their nearby users, thus offloading traffic from 
often congested cellular networks [3]-[6]. This is facilitated 
by the spatiotemporal correlation in the content demanded i.e., 
repeated requests for the same content from different users 
across various time instants [7]-[9]. Storing the popular files 
at the “network edge”, such as in small cells, switching centers 
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or handheld devices, termed caching, offers an excellent way 
to exploit this correlation in the content requested by the 
users [10]-[15]. Cache-enabled D2D networks are attractive 
due to the possible linear increase of capacity with the number 
of devices that can locally cache data [4], [16]-[21]. 

The performance of cache-enabled D2D networks funda¬ 
mentally depends upon i) the locations of the devices, and ii) 
how is cache placed on these devices. For instance, consider 
the device-centric placement where the content is placed on a 
device close to the particular device that needs it. While this 
is certainly beneficial for the device with respect to which the 
content is placed, high performance degradation can happen 
if another device in the network wants to access the same 
content from the device on which it was cached. As a result, 
we focus on the cluster-centric placement, where the goal is 
to improve the collective performance of all the devices in the 
network measured in terms of the coverage probability and 
area spectral efficiency (ASE). This requires several new results 
for the coverage probability and ASE, where the receiving 
devices of interest and/or the devices that contain content of 
interest for these receivers are parametrized in terms of the 
cluster-center. These results are the main focus of this paper. 

A. Related Work and Motivation 

Existing works on the modeling and analysis of D2D 
networks has taken two main directions. The first line of 
work focuses on characterizing the asymptotic scaling laws for 
cache-enabled D2D networks using the well-known protocol 
model, see [17]-[19], [21] for a small subset. These works rely 
on a grid-based clustering model where the space is tessellated 
into square cells with devices in each cell forming distinct 
clusters. While these works provide several key insights, the 
key limitation is the use of the protocol model, which assumes 
that the communication between two nodes is possible only if 
the intended receiver is: i) within collaboration distance of the 
intended transmitter, and ii) outside the interference range of 
all other simultaneously active transmitters [22]. The second 
line of work considers the so-called physical model, where 
the successful communication between two nodes is based 
on the received signal-to-interference-and-noise ratio, unlike 
the protocol model where it is simply based on the distance 
[22]. Tools from stochastic geometry have been used for the 
tractable characterization of the key physical layer metrics, 
such as the rate and coverage [23]-[25]. These tools have 
resulted in significant advancement in the tractable modeling 
and analysis of downlink and uplink cellular networks [26]- 
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[29]. Motivated by this encouraging progress, there has re¬ 
cently been a surge of interest in applying these tools to the 
analysis of D2D networks. These works are discussed next. 

Depending upon the spectrum allocated for D2D transmis¬ 
sions, the D2D networks can be classified into two categories: 
in-band and out-of-band. In the in-band D2D networks, D2D 
and cellular networks coexist in the same spectrum. Using 
tools from stochastic geometry, several coexistence aspects, 
such as mode selection between cellular and D2D [30]- 
[32], coexistence of D2D and unmanned aerial vehicle [33], 
distributed caching in D2D networks [34], and D2D interfer¬ 
ence management to protect cellular transmissions [35]-[39], 
have been addressed. On the other hand, in the out-of-band 
D2D, as the name implies, orthogonal spectrum is allocated 
for D2D and cellular transmissions. For this setup, various 
aspects and applications of D2D networks, such as multicast 
D2D transmissions [40], D2D communication with network 
coding [41], and traffic offloading from cellular networks to 
D2D networks [42], have been studied. 

To lend tractability, the common approach in all the above 
mentioned stochastic geometry-based works for D2D networks 
is to model the locations of D2D transmitters (D2D-Txs) as 
a Poisson Point Process (PPP), and the locations of D2D 
receivers (D2D-Rxs) via two approaches: i) D2D-Rxs lie at a 
fixed distance from their intended D2D-Txs [30], [33], [35]- 
[39], and ii) D2D-Rxs are uniformly distributed in a circular 
region around their intended D2D-Txs [31], [41], [42]. While 
these models provide several useful design insights, they suffer 
from a key shortcoming of not being able to capture the notion 
of device clustering, which is quite fundamental to the D2D 
network architecture [17]-[19], [43]-[45]. This shortcoming 
was addressed in our very recent work [46], [47], where we 
modeled the device locations by a Poisson cluster process to 
analyze the performance of device-centric content placement 
policies. In contrast, the current work takes a cluster-centric 
approach where the focus is on placing content so as to 
optimize the performance of the whole cluster rather than the 
individual devices. The analysis involves the characterization 
of several new distance distributions in Poisson cluster pro¬ 
cesses. More details along with other main contributions are 
explained in detail below. 

B. Contributions and Outcomes 

Tractable model for cache-enabled D2D networks: We 
develop a realistic analytic framework to study the perfor¬ 
mance of cluster-centric content placement policies in a cache- 
enabled D2D network. Modeling the locations of the devices 
by a Poisson cluster process (in particular a slight variation of 
a Thomas cluster process) and using tools from stochastic ge¬ 
ometry and stochastic orders, we first prove that the collective 
performance of all the devices in a given cluster in terms of 
coverage probability is improved when the content of interest 
for each device is placed at the device closest to the cluster 
center. This policy, however, may not always be feasible due to 
the limited storage capacity and/or the energy of the closest 
device. Besides, placing all the content on a single device 
limits frequency reuse within a cluster to one, which may 


Table I: Summary of notation 


Notation 

Description 

<E>c;Ac 

An independent PPP modeling the locations of D2D cluster center, 
density of D2D cluster center 

X 

The location of cluster center 

a, b 

The relative location of cluster member form cluster center 

R, r 

The serving distance, where a realization of R is denoted by r 


Set of devices inside the cluster 


Set of possible transmitting and receiving devices 

N,Nt,Nr 

Number of total, possible transmitting and receiving devices 

A^] rfia,, B^-, mb 

Set of simultaneously active devices inside the cluster with mean ma and mb 


Scattering variance of cluster member location around cluster center 

Pi 

Transmit power of devices engage in D2D communications 

Q. 

Path loss exponent corresponding to the D2D link; a > 2 

hax 1 

Exponential fading coefficients with mean unity 

P 

SIR threshold for successful demodulation and decoding 

Pc 

Coverage probability 

ASE 

Area spectral efficiency 


not be optimal in terms of the network throughput. As a 
result, we explore more general scenarios in which the content 
is distributed across devices in the cluster. The analysis of 
such scenarios require new methodology where the receiving 
devices of interest and/or the devices that contain their content 
of interest are parameterized in terms of their location relative 
to the cluster center. This forms the main technical contribution 
of the paper. More details are provided next. 

ASE and coverage probability analysis: We derive easy 
to use expressions for coverage probability and ASE for the 
following three general cases: (i) k-Tx case: receiver of interest 
is chosen uniformly at random in a cluster and its content of 
interest is available at the closest device to the cluster 
center, (ii) f-Rx case: receiver of interest is the closest 
device to the cluster center and its content of interest is 
available at a device chosen uniformly at random from the 
same cluster, and (iii) baseline case: the receiver of interest 
is chosen uniformly at random in a cluster and its content 
of interest is available at a device chosen independently and 
uniformly at random from the same cluster. A common feature 
in these cases is the parameterization of the receiver of interest 
and/or its serving device in terms of its location relative to 
the cluster center. These results provide key insights into the 
performance of cluster-centric content placement policies. A 
key intermediate step in the analysis is the characterization of 
distances from the D2D-Rx of interest to its serving device, 
and intra- and inter-cluster interfering devices for these cases. 

New generative model for cluster-centric D2D networks: 
The analysis of fc-Tx and f-Rx cases described above shows 
that the network performance improves significantly when the 
device(s) on which the content is cached or the device(s) 
requesting content from cache are biased to lie closer to the 
cluster center. This means that besides the D2D link of interest, 
the intra-cluster interfering links may be more likely to have 
a transmitter or receiver closer to the cluster center. To study 
the effect of this behavior on the network performance, we 
propose a generative model in which the device locations 
follow a double-variance Thomas cluster process, where each 
cluster consists of a denser and a sparser subcluster. Sampling 
the locations of the transmitters or receivers uniformly at 
random from the denser subcluster allows us to model the 
above described biasing behavior fairly accurately. 
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II. System Model 

We consider a clustered D2D network where the content 
of interest for devices of a given cluster is cached in the 
same cluster. This is inspired by the fact that the popular 
content may vary significantly across clusters. For instance, 
users in a library may be interested in an entirely different 
set of files than the users in a sports bar. Besides, larger 
inter-cluster distances make it difficult to establish direct 
communication across clusters. Note that while our model is, 
in principle, extendible to include inter-cluster communication, 
we will limit our discussion to more relevant case where 
direct communication is only between two devices of the 
same cluster. More details on how the content is placed in 
the devices of a given cluster will be provided in Section III. 
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Figure 1: Illustration of D2D cluster network when cluster members 
(devices) are normally distributed around cluster center with = 40. 


A. System Setup and Key Assumptions 


We model the locations of the devices by a Poisson cluster 
process in which the parent points are drawn from a PPP <I>c 
with density Ac and the offspring points are independent and 
identically distributed (i.i.d.) around each parent point [48]. 
The parent points and offsprings will be henceforth referred 
to as cluster centers and cluster members (or simply devices), 
respectively. The cluster members (or devices) around each 
cluster center a; € $c are sampled from an i.i.d. symmetric 
normal distribution with variance in Therefore, the 
density function of the location of a cluster member relative 
to the location of its cluster center, a € is 


/a( a) = 


1 


2T:al 


exp - 


2a2 


( 1 ) 


If the number of cluster members in each cluster is Poisson 
distributed, this setup corresponds to the well-known Thomas 
cluster process [49]. Note that we will put some restrictions on 
the number of cluster members to facilitate characterization of 
distance distributions in the sequel. Therefore, our setup can 
be interpreted as a variant of Thomas cluster process. 

Denote the set of devices belonging to the cluster centered 
as X S 4)c by Partition this set into two subsets of (i) 
possible transmitting devices denoted by A/]f, and (ii) possible 
receiving devices denoted by Afff Within each cluster, the set 
of simultaneously active transmitters is denoted by 
and hence the set of simultaneously active transmitters in the 
whole network can be expressed as: 




To keep the model general, we assume that the number of 
simultaneously active transmitters \A^\ is not necessarily the 
same for each cluster. More specifically, |Al“| is modeled as a 
Poisson distributed random variable with mean rria. 

Without loss of generality, we focus on a randomly chosen 
cluster, termed representative cluster, with its cluster center 
denoted by Xq S $(,. For this cluster, we assume that the 
total number of devices is = N, and the number of 

possible transmitting devices is = Nt- This assumption 

is made to facilitate order statistics arguments that appear 
in the characterization of distance distributions in the sequel. 
Note that for a meaningful analysis, the link corresponding 


to the D2D-Rx of interest in the representative cluster needs 
to be active. Once the location of the D2D-Tx of interest 
is fixed, the set of other simultaneously active transmitters 
in the representative cluster is sampled uniformly at random 
from the remaining iVt — 1 positions. Therefore, it is assumed 
that the number of intra-cluster interfering devices is Poisson 
distributed with mean fh-a — 1 conditioned on the total being 
less than A^t — 1- As a result, the average number of active 
devices in the representative cluster is rfia, which is consistent 
with the assumption made above regarding the number of 
simultaneously active transmitters per cluster. 


B. Channel Model 

Recall that the cluster center of the representative cluster 
is assumed to be located at xg G and hence the D2D-Rx 
of interest belongs to the set . Without loss of generality, 
the analysis is performed at the D2D-Rx of interest in the 
representative cluster, which is assumed to be located at the 
origin. The D2D-Txs are assumed to transmit at the constant 
power Pd- The content of interest for the D2D-Rx of interest 
is available at the device located at oo+xq, where og indicates 
the location of this device relative to cluster center xg. Hence, 
the received power at the D2D-Rx of interest is 


P = PdhaoWxo + aoW 


where hag ~ exp(l) models Rayleigh fading, and a is the 
power-law path loss exponent. The total interference experi¬ 
enced by the D2D-Rx of interest can be written as a sum 
of two independent terms. First, the interference from the set 
of devices inside the representative cluster, say intra-cluster 
interference, is given by 

X! -Pd^a^J|a:o + a||““. (2) 

aGA“=o\ao 

Second, the interference from the devices outside the repre¬ 
sentative cluster, say inter-cluster interference, is given by 

/inter = J|a: + a||(3) 
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Now, the signal-to-interference ratio (SIR) at the D2D-Rx of 
interest at a distance i? = ||a:o + ao|| from the serving device, 
where a realization of R is denoted by r, is: 

Pi h r ^ 

SIR(r) = T. -• 

Winter “r -tintra 

Note that the SIR expression is not a function of the transmit 
power Pd and therefore without loss of generality, we assume 
that Pd = 1. For this setup, we study network performance 
in terms of coverage probability and ASE which are formally 
defined next. 

Definition 1 (Coverage probability). The probability that SIR 
of an arbitrary link of interest at the receiver exceeds the 
required threshold for successful demodulation and decoding. 

Pc =E[l{SIR(r) >/3}], (5) 

where (3 is a pre-determined threshold for successful demod¬ 
ulation and decoding at the receiver. 

Definition 2 (Area spectral efficiency). The average number 
of bits transmitted per unit time per unit bandwidth per unit 
area can be defined as: 

ASE = Alog 2 (l + /3)E[l{SIR(r) > /3}], (6) 

where A is the number of simultaneously active transmitter 
per unit area. 

III. Coverage Probability and ASE 

This is the first main technical section of the paper, where 
we first characterize the coverage-optimal cluster-centric con¬ 
tent placement policy. We show that under this policy the 
content of interest for all the devices must be stored at the 
device closest to the cluster center. However, this may not be 
feasible due to storage and/or energy constraints of mobile 
devices. Besides, such a placement would limit the number of 
simultaneously active D2D connections over a given frequency 
band in a given cluster to one. Since aggressive frequency 
reuse is one of the advantages of D2D, this is clearly not 
desirable. As a result, we assume that the content is distributed 
across devices in a cluster. To enable the analysis of cluster¬ 
centric content placement policies for this setup, we define 
several cases of fundamental interest where the D2D-Rx of 
interest and/or the device which has its content of interest are 
parameterized in terms of their locations with respect to the 
cluster center. Easy-to-use expressions for coverage probability 
and ASE are then derived for these cases. 

A. Coverage-Optimal Content Placement 

In this subsection, we study the coverage-optimal content 
placement problem in the proposed clustered D2D model. 
Note that while it is preferable to place the content required 
by each device at its immediately neighboring device, such 
device-centric content placement is not realistic. Therefore, 
we focus on the cluster-centric content placement, where the 
goal is to place the content in such a way that it improves 
the collective performance of the whole network. This can be 


achieved by fixing the point of reference for content placement 
to be the cluster-center instead of a particular receiver. To 
fix this key idea, we begin with a simple problem where we 
assume that the content of interest for the whole cluster is 
placed at a single device in Aff° that is closest to the 
cluster center, where k = 1 and k = Nt correspond to the 
closest and farthest devices from the set to the cluster 
center, respectively. Our first goal is to find the value of k 
that optimizes the performance of the whole cluster. We cast 
this problem as the coverage maximization problem, where 
coverage probability of a D2D-Rx of interest is 

Pc=E[l{SIR(||xo + Sfc||) >/3}], (7) 

where xq is the location of the cluster center, and Sk is the 
location of the closest device to the cluster center, which is 
also the serving device. The optimal value of k that maximizes 
this coverage probability is derived in the next Lemma. 

Lemma 1. The optimal value ofk that maximizes the coverage 
probability given by (7) for the whole cluster is 

arg 2 ....,At} E[1 {SIR(||xo -f Sfc||) > /3}] = 1. (8) 

Proof. See Appendix A. 

■ 

An intuitive interpretation of the above result is that all the 
devices in a given cluster should be served by a device that 
is on an average closest to all of them. As proved formally 
in the above Lemma, this device is the one that is closest 
to the cluster center. While this result is potentially useful in 
determining the coverage-optimal location of a cache-enabled 
small cell or a dedicated storage device, this simple policy 
limits the frequency reuse capability of D2D networks by 
concentrating all the content at a single device. Besides, such a 
policy may be infeasible due to storage and energy constraints 
of mobile devices. Therefore, it is important to distribute the 
content across multiple devices in a cluster. 

As noted above, for cluster-centric content placement, the 
point of reference will be the cluster-center instead of a 
particular receiver. This means the D2D-Rx of interest and/or 
the device that has cached its content of interest can be 
parametrized in terms of their locations with respect to the 
cluster center. Lor instance, it is precise to say that a receiver 
of interest will have its content of interest cached at a device 
that is closest to the cluster center from the set 
where k G The value of k will depend upon the 

content placement strategy being adopted, as discussed in the 
context of optimizing the total hit probability in Section V-B. 
Lor performance comparison, we also consider a random 
placement strategy, where the content requested by the D2D- 
Rx of interest is available at a device chosen uniformly at 
random from the set Aff°. 

To analyze the performance of the above setup, we also need 
to define how the D2D-Rx of interest is chosen. We consider 
two choices: (i) D2D-Rx of interest is parameterized with 
respect to the cluster center as done for the transmitter above 
(say closest to the cluster center from the set and (ii) 

the D2D-Rx of interest is chosen uniformly at random from 
the cluster. While the latter provides insights into the typical 
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network performance, the former is useful in understanding 
how the performance of devices located towards the center 
of the cluster (small values of i) differs from those located 
at the edge of the cluster (large values of £). For this setup, 
we focus on the following three cases, each providing useful 
insights into the performance of D2D networks: 

• k-Tx case: In this case, the D2D-Rx of interest is chosen 
uniformly at random from a given cluster and the content 
of interest for this receiver is available at the closest 
transmitting device to the cluster center (in the set 
from the same cluster. By tuning the value of k, we 
can study the effect of the location of the content/cache 
(relative to the cluster center) on the performance. 

• £-Rx case: In this case, the D2D-Rx of interest is the 

closest device to the cluster center in the set 
and its content of interest is available at a device chosen 
uniformly at random in . By tuning the value of f, 
we can understand how the performance of users located 
towards the center of the cluster differ from those located 
towards the cluster edge. 

• Baseline case: In the baseline case, we assume that the 
D2D-Rx of interest is chosen uniformly at random from 

and the device containing its content of interest is 
also chosen uniformly at random from . This simple 
case will act as a baseline for performance comparisons. 
Note that we can, in principle, define fc-Tx f-Rx case, where 
the D2D-Rx of interest is the closest device to the cluster 
center from the set and its content of interest is available 
at the k^^ closest device to the cluster center from the set . 
Due to lack of space and the fact that the essence of this case 
will be captured approximately in the new generative model 
studied in Section IV, we do not consider this explicitly. 

As discussed in detail in Section II, the effect of frequency 
reuse is studied by assuming that multiple D2D links in a 
given cluster can be activated simultaneously. In particular, 
once the location of the serving device is decided, the locations 
of intra-cluster interfering transmitters are sampled uniformly 
at random from the remaining points of . Note that the 
location of these interfering devices can be sampled in more 
sophisticated ways (e.g., biased to lie closer to the cluster 
center). This will be discussed in Section IV. We now derive 
the coverage probabilities for the three cases described above 
in the next subsection. 

B. Coverage Probability Analysis 

Before going into the detailed analysis of coverage proba¬ 
bility, we characterize the distributions of the distances from 
intra- and inter-cluster devices to the D2D-Rx of interest under 
various polices. Using these distance distributions, we derive 
the Laplace transform of distribution of intra- and inter-cluster 
interference distributions. As will be evident from our analysis, 
characterizing Laplace transform of interference distribution is 
key intermediate result for the coverage probability analysis. 
The distances between the D2D receiver of interest and the 
various inter/intra cluster interfering devices are in general 
correlated. Focusing first on the intra-cluster devices, denote 
the distances from the D2D-Rx of interest to the intra-cluster 


devices by {rw}, where w = ||a;o + a|j. Clearly these distances 
are correlated because of the common distance between the 
cluster center to the D2D-Rx of interest, vq = ||a:o||- As 
discussed in detail in [46] for device-centric content place¬ 
ment, this correlation can be handled by conditioning on the 
common distance r^o = lla^oll^ after which the distances {re} 
become conditionally i.i.d., which lends tractability to the 
analysis of the Laplace transform of interference distribution, 
thus resulting in tractable expressions for coverage probability 
and ASE. For the current cluster model, the conditional distance 
distribution fw{w\i:o) is characterized by Rician distribution 
[46]: 


fw{w\iyo) = -ly exp 


-f t'o \ 
) 



w > 0. 

(9) 


where Iq{.) is the modified Bessel function of the first kind 
with order zero. Since probability density function (PDF) of 
the Rician distribution will be frequently used in the sequel, 
we define its functional form below to simplify the notation. 


Definition 3. (Rician distribution). The PDF of the Rician 
distribution /y(y|^) is 


Ricepdf(y,z;a2) = ^ exp Jp (p) , 

where a is the scale parameter of the distribution. 


2 / > 0 , 
( 10 ) 


Note that while the intra-cluster distances to the D2D- 
Rx of interest become conditionally i.i.d. by conditioning on 
^0 = II3^0II> there is still a possibility of dependence between 
the distances to the serving and interfering devices when the 
serving device is not chosen uniformly at random from the 
cluster. The easiest way to understand this dependence is 
by recalling that while the sequence of distances to intra¬ 
cluster devices, w = ||xo -f a||, is conditionally i.i.d., the 
“ordered” choice of the serving device impacts the distribution 
of the remaining elements in the sequence (distances to the 
interfering devices). This is clearly true in the /c-Tx case, 
where the serving device is chosen to be the k*^ closest device 
to the cluster center. However, it turns out that this dependence 
can be handled by first conditioning on the distance from the 
serving device to the D2D-Rx of interest, denoted by tk = 
||sfe||, and then partitioning the intra-cluster interfering devices 
into two subsests: (i) devices that are closer than the serving 
device to the cluster center, denoted by a € where the 
distance to the cluster center is denoted by tin = II^ll^ (ii) 
devices that are farther than the serving device to the cluster 
center, denoted by a € where the distance to the cluster 

center is denoted by tout = ||o||- Please refer to Fig. 2 for the 
pictorial representation. In the following Lemma, we prove 
that the distances from devices in a € A^u a € Aout 
respectively i.i.d., which lends tractability to the interference 
analysis. This result along with the conditional distribution of 
the distances is given next. 

Lemma 2 (Distance of intra-cluster interfering device to the 
D2D-Rx of interest in the fc-Tx case). The distances from the 
intra-cluster interfering devices to the D2D-Rx of interest in 
the k-Tx case, i.e., {w = ||a:o + a||} are conditionally i.i.d.. 
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Figure 2: Illustration of intra-cluster devices for the fc-Tx case. 


conditioned on vq = ll^ioll. ond t = ||a|| (where t can be either 
fin or foutj. with PDF 


1 


fw{w\vo,t) = - 


where, if t = f, 


w/uot 




r, \l^ 0 -t\ <W < VQ+t, 


7 i'in ^ i'k 


/Ti„(fin|ffc) = i i_exp(^-^^ 

^in ^ 

if t = tout, then. 


/T„„t (t out \tk) = 


7 ^2 \ ? ^out > tk 


0 , 


fout ^ ffc 


( 11 ) 


( 12 ) 


(13) 


Proof. See Appendix B. 


C — 1 (rtta — l)^e r/s s I \ ■ 

i = hj=0 - i!-’ i(l +gm) IS 

regularized incomplete beta function, and density functions 
of fwiwWo, t), fTiWfn\tk), and /To„t (*out|ffc) are given by 
Lemma 2. 

Proof See Appendix C. ■ 


For the other two cases (£-Rx and baseline), the selection 
of serving device is done uniformly at random, which does 
not induce any dependence in the distances from the serving 
and interfering devices, leading to a simpler expression for 
the Laplace transform of intra-cluster interference distribution 
for these cases. Note that conditioning on uq = ||a;o|j is still 
necessary to handle the correlation induced by the common 
term xo in the intra-cluster distances, as discussed earlier in 
this section. 


Lemma 4. For the £-Rx and baseline cases, the conditional 
Laplace transform of the intra-cluster interference distribu¬ 
tion, conditioned on the distance from the D2D-Rx of interest 
to the cluster center, vq = ll^ioll. is 


Nt-l 


WsWo)= (m{wWo)) 




n—0 


2^ng-(ma-l) 

n\^ 

(15) 


with M(wWo) = /o°° i+L— fv/{wWWjA.w. Assuming < 
Ni, we have 

^/i„tra(sko) ^exp - (rfia - 1)(1 - M(^wWo))), (16) 

where ^ ^ ° fw{wWo) = 

Ricepdf (w, uq] a^). 

Proof See Appendix D. ■ 


Using this distance distribution, the exact expression of 
the conditional Laplace transform of intra-cluster interference 
distribution in the fc-Tx case is given next. Please note that the 
corresponding result appearing in the shorter version of this 
paper [1] is an approximation. 


Lemma 3. In the k-Tx case, the conditional Laplace transform 
of distribution of intra-cluster interference (2), conditioned on 
uq, where content of interest is placed at distance tk = ||sfc|| 
from cluster center is = 


Nt-l Sm 

EE 

n =0 i =0 


(?y(i-p)' 


-I 


/(l-p;n-5m,l+5m) 


Minis,tkWW'' 


^ Miout {.S, tk\ Uq ) 


nW 


with. 


(14) 


Remark 1. As discussed in the sequel, the optimal number of 
simultaneously active links is much smaller than the total num¬ 
ber of potential transmitters. Hence, the assumption rfin <§; 
taken to derive the simpler expression (16) in Lemma 4 is fairly 
reasonable. As a result, (16) will in general be quite accurate. 

Remark 2. While Lemma 4 is exact for i-Rx and baseline 
cases, it also provides a tractable approximation for the k- 
Tx case, whose exact expression for the Laplace transform 
of intra-cluster interference distribution given by Lemma 3 is 
much more complicated due to the presence of two summa¬ 
tions. Lemma 4 is an approximation for the k-Tx case because 
it ignores the effect of "ordered” selection of serving device 
on the distance distributions of the intra-cluster interference 
devices. However, this approximation for the k-Tx case will 
be numerically shown to be quite tight in Section V. 


Mi„(s,ffc|r'o) = [ [ " E Wr, Iffc)dwdfin, 

Jo J™p„ 1 + 

■^^out (s, ifc ^o) - If T . ^ ^Tout (^out 

where = W - fn\, = vo + Un, = Wo - iout|, 

Wont = + tout, P = Ym’ 5 m = mm{n,k - 1), 


The Laplace transform of intra-cluster interference distribu¬ 
tion given by Lemma 4 can be simplified further under the 
following assumption without loosing much accuracy. 

Assumption 1 (Un-correlated intra-cluster distances assump¬ 
tion for fc-Tx and baseline cases). Since the devices are 
normally distributed around the cluster centers, the distances 
from the intra-cluster devices to the D2D-Rx of interest are 
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Rayleigh distributed with the following PDF when the D2D- 
Rx of interest is chosen uniformly at random [46] 


w 

/»■(.») = 5 ^ exp — 


w > 0. 


(17) 


However, as discussed earlier in this section, the distances 
are correlated due to the presence of the common distance 
1^0 = ||a;o||. to which we conditioned on this distance in 
(9). However, if we ignore this correlation, we can simplify 
the analysis by assuming that the distances are i.i.d. Rayleigh 
distributed with the PDF given by (17). This approximation is 
however not applicable for the £-Rx case where the D2D-Rx 
of interest is not chosen uniformly at random. 

Under this assumption the approximation for Laplace trans¬ 
form of intra-cluster interference distribution is given next. It 
is applicable for the fc-Tx and baseline cases. 

Corollary 1. Under Assumption 1, the Laplace transform 
of intra-cluster interference distribution in k-Tx and baseline 
cases is 


'C/i„tra(s) = ^ exp ( - (ma - 1) 




sw 


1 + sw~^ 


fwiw)dv?j, (18) 


where ^ = J2f=o —" ond fw{w) given by (17). 

We will use this simpler expression to provide easy to 
compute expression for coverage probability later in this 
section. We now derive the Laplace transform of inter-cluster 
interference distribution. Recall that the inter-cluster interferers 
are sampled uniformly at random in all three cases, which 
means the following result is exact for all three cases. 

Lemma 5. For all three cases, the Laplace transform of 
distribution of inter-cluster interference at D2D-Rx of interest 
in (3) is 


L, 


Xs) = exp (-2TrXc J (l-exp (-rrii, J 


su 
1 + 


fu{u\iy)du^udu^y (19) 


where fu{'<Av) = Ricepdf (u, v] a^). 

Proof The Laplace transform of the inter-cluster interference 
distribution is special case of the Lemma 8. Since the proof 
follows on the same lines as that of Lemma 8, its skipped. ■ 


1) Coverage probability analysis of k-Tx case: Recall that 
the D2D-Rx of interest in this case is chosen uniformly 
at random and the D2D-Tx of interest is the closest 
transmitting device to the cluster center (in the set We 

first derive the serving distance distribution for this case. 


Lemma 6. The PDF of the serving distance, i.e., r = ||a;o + 
Sfc||, conditioned on the distances vq = ||xo|| and tk = ||sfc|| 
for the k-Tx case is 

fRirWo,tk) = - I - - t^\ <c r < uq-\- tk, 


( 20 ) 

with 

fvofyo) = ^ exp . *70 > 0 (21) 

= (fc - 1 )k1 - ~ F{tk)f^-'^f{tk) 

( 22 ) 

where flfk) = ^ exp(-^), and F{tk) = 1 - exp(-^). 

Proof The PDF of serving distance r = ||a:o+Sfc|| conditioned 
on the Vo and tk, i.e., fR{r\iyo,tk) can be derived exactly 
on the same lines as fwiw\Ro,t) given by (11). Hence, the 
proof is skipped. Here, fvg (r'o) is Rayleigh distributed owing 
to the fact that the D2D-Rx of interest is a randomly chosen 
device where devices are normally scattered around the cluster 
center. Finally, for note that the distances of intra-cluster 
devices to the cluster center are i.i.d. Rayleigh distributed with 
Tk being the smallest sample out of Nt elements, whose 
distribution follows by order statistics (see [50, eq (3)]). ■ 

Using this result, we now derive the coverage probability 
for the k-Tx case in the following theorem. 

Theorem 1 (Coverage probability: fc-Tx case). Using Laplace 
transform of distribution of interference in (14), and (19), the 
coverage probability of the D2D-Rx of interest is 

^oo ^oo 

t/ 0 ^0 J 

X fvoMfTk{tk)drdiyodtk, (23) 

with = |u — tk\, and = v -\- tk, where fR{r\vo,tk), 
/v'o(^o). ond /Tfc(ffc) are given by (20), (21), and (22) respec¬ 
tively. 


Proof. From the definition of coverage probability, we have 


ETfcEvoEi? 

P| 

^ (-Winter ^intra) 

R,Vo,Tk}] 

ErfcEyoEfl 

E 

exp ( /3r {-/^inter -^intra)) 

R,Vo,Tk]] 


where (a) follows from hoxo ~ exp(l). The result follows 
from the fact that intra- and inter-cluster interference powers 
are independent, followed by the expectation over R given vo 
and tk, followed by expectation over Vq and Tk. The PDFs of 
Vo and Tk are given by (21) and (22), respectively. ■ 

As discussed in Remark 2, the exact expression for the 
Laplace transform of the intra-cluster interference distribution 
given by (14) in Lemma 3 is quite complicated due to the 
presence of two summations. To improve tractability, the 
simpler expression of Lemma 4 can be used. This leads to an 
approximation since the dependence of the distances from the 
intra-cluster interfering devices on the selection of the serving 
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device is not captured. The approximate result is given next. 
The proof follows on the same line as that of Theorem 1. 

Corollary 2. Using Laplace transform of intra-cluster inter¬ 
ference distribution given by Lemma 4, the coverage proba¬ 
bility of k-Tx case can be approximated as 

poo poo poo 

Jo Jo Jo 

fvoMfn{tk)drdvodtk, (24) 

where is given by (19), and /i?(r|r'o, 4), fvoi^o), 

fTkitk) are given by (20), (21), and (22) respectively. 


at xq + a where xg = si, and a is sampled from zero-mean 
complex Gaussian random variable. ■ 

3) Coverage probability analysis of the baseline case: In 
the baseline case, we assume that both the D2D-Rx of interest 
and D2D-Tx of interest are sampled uniformly at random. 
The coverage probability for this case is given next. For the 
complete proof, please refer to [46], where the same case was 
used as the baseline case for device-centric content placement 
strategies. Here we just provide a proof sketch. 

Theorem 3 (Coverage probability: baseline case). The cover¬ 
age probability in the baseline case is 


Although the above coverage probability expression for k- 
Tx case seems to be involved, it can be easily evaluated by 
Quasi-Monte Carlo numerical integration methods (because 
the integrations are essentially expectations) [51]. Using the 
approximation of the Laplace transform of the intra-cluster 
interference distribution given by Corollary 1, we can further 
simplify coverage probability expression in the next corollary. 


Corollary 3. By ignoring intra-cluster distance correlations 
under Assumption 1, the coverage probability of the k-Tx case 
can be approximated as 


pOO /*00 


pTx ^ 

Cfc — 


fo Jo 


(25) 

where fR{r\tk) = Ricepdf (r, 4; ct^), and fr^itk) is given by 

( 22 ). 


P? = 


^h^t„{dr°‘)^h„trSdr°‘Wo)fR(r\vo)fvo(vo)drAvQ, 

( 28 ) 


where fR{r\vg) = Ricepdf (r, r-o; 0 -^), and /vo(t'o) = 
^exp (-^). 

Proof The proof follows on the same lines as Theorem 1. By 
definition, the coverage probability is 




E 


exp ( /3r (/inter “t” /intra)) 


R,Vg 


(29) 


where the PDF of serving distance r = ||a:o + ao|| conditioned 
on ug = ||xo|| is Rician distributed [46]. Further, the D2D-Rx 
of interest is chosen uniformly at random, which is sampled 
from a Gaussian distribution in and hence ug = ||a;o|| 
simply follows Rayleigh distribution. ■ 


Proof See Appendix E. 


The tightness of the approximation will be validated in the 
numerical results section (Section V). 

2) Coverage probability analysis of £-Rx case: We now 
derive the coverage probability for the /-Rx case, where the 
D2D-Rx of interest is the closest device to the cluster 
center from the set Aff° and its serving device is chosen 
uniformly at random from the set . 

Theorem 2 (Coverage probability: f-Rx case). The coverage 
probability in £-Rx case is 


P 


Rx 



poo 

/ iPr‘^\L)fR{r\te)fT, {te)drdte, 

Jo 

( 26 ) 


with friite) = 

where f{te) = ||exp(-^), F(Q) = 1 - exp(-^), and 
/it(r|4) = Ricepdf(r,4;cra). 


Proof The proof follows on the same lines as the proof of 
Theorem 1. By definition, the coverage probability is 


pRx 

Ci 


Et,E_r 


E 


exp ( fiv (/inter 4 /intra)) 


R, Ti 


(27) 


where the PDF of serving distance r = ||s^ -f a|| conditioned 
on ti = 11 Sr 11 can be derived on the same lines as the PDF 
of the serving distance in Corollary 3. This is because in this 
case, file of interest is available inside the cluster uniformly at 
random and D2D-Rx of interest is closest receiving device 
to the cluster center. Thus, the D2D-Tx of interest is located 


C. Area Spectral Efficiency Analysis 

After studying the coverage probability for all three policies 
in the previous subsection, we now focus on the area spectral 
efficiency (ASE), which is defined in Definition 2. This defini¬ 
tion is specialized to our setup in the following proposition. 

Proposition 1. The ASE/or the three cases is given by: 

ASE = maAclog2(l+/3)Pc, (30) 

where Pc is given by (23), (26), and (28) for k-Tx, £-Rx, and 
baseline cases respectively. Here, rhg^Xc denotes the average 
number of simultaneously active transmitting devices. 

Remark 3 (Trade-off between channel orthogonalization and 
more aggressive spectrum reuse). Intra-cell channel orthog¬ 
onalization, i.e., a small number of simultaneously active 
devices per cluster, reduces intra-cluster interference at the 
expense of less spectrum reuse. We cast this classical tradeoff 
between higher number of simultaneously active links (i.e., 
more spectrum reuse) and higher interference as the problem 
of finding the optimum rfig. that maximizes ASE: 

ASE* = max WaAc log 2 (l +/3)Pc- (31) 

We will revisit this trade-off over the number of simultaneously 
active links in the numerical results section. By solving this 
ASE optimization problem numerically, we will demonstrate 
that optimum number of simultaneously active links is signif¬ 
icantly different for the three cases, which means it is highly 
dependent on the choice of content placement policy. 
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Figure 3: Illustration of the double-variance Thomas cluster model 
when cluster members (devices) are normally distributed around the 
cluster center with Ca = 10 and crb = 30. 

IV. New generative model eor the cluster-centric 
D2D NETWORKS 

A key takeaway from the analyses of fc-Tx and £-Rx cases, 
which will become more apparent in the numerical results 
discussed in Section V, is that the network performance 
improves signihcantly when the device(s) on which the content 
is cached or the device(s) requesting content from the cache 
are biased to lie closer to the cluster center. This means that in 
addition to the D2D link of interest, the intra-cluster interfering 
links may be more likely to have a transmitter or receiver 
closer to the cluster center. Incorporating such a behavior in 
the original model of Section II will require hxing the indices 
of the interfering devices in each cluster (relative to the cluster 
centers), as we did for the serving and receiving devices in 
the fc-Tx and f-Rx cases. While this is certainly doable in 
principle, the final expressions will be prohibitively complex 
due to the dependence amongst all the distances involved in 
the coverage probability evaluation. For instance, the intra¬ 
cluster distances will have to be jointly handled through their 
joint distribution that will be evaluated using order statistics 
on the same lines as the serving distances in fc-Tx and ^-Rx 
cases. Deconditioning on such joint distributions will result in 
multi-fold integrals that are not easy to evaluate. 

Therefore, to study the effect of this biasing of potential 
transmitters and receivers towards the cluster center, we pro¬ 
pose a generative model in which the device locations follow 
a double-variance Thomas cluster process, where each cluster 
consists of a denser and a sparser subcluster. As discussed in 
this section, selecting the locations of the transmitters or re¬ 
ceivers uniformly at random from the denser subcluster allows 
us to model the above described biasing while maintaining 
tractability. The generative model is illustrated in Fig. 3. More 
formal details about the proposed model are presented next. 

A. System Setup and Channel Model 

We model the clustered D2D network using a double¬ 
variance Thomas cluster process whose cluster centers are 
distributed according to a PPP $c of density Ac. A cluster 
around x S $c is formed of two subclusters, denser and 


sparser, of normally distributed devices with scattering vari¬ 
ances cr^ and (7^’ respectively. The analysis will be performed 
at a typical device, i.e., a device chosen uniformly at random 
from the either subclusters. This means that the D2D-Rx of 
interest will be normally distributed around the cluster center 
with variance or a^. The simultaneously active transmitters 
of the two subclusters are denoted by (denser) and 
(sparser), where \A^\ and \B^\ are Poisson distributed with 
mean Wa and mb, respectively. Since we want to bias the 
location of the D2D-Tx of interest towards the cluster center, 
we sample it uniformly at random from the denser subcluster 
A^°. While the other case in which it is sampled from B^° 
is not important for the current discussion, it can be handled 
exactly on the same lines. As was the case in the original 
model, the number of intra-cluster interfering devices in A^° is 
modeled as a PPP with mean Wa — 1 to ensure that the average 
number of simultaneously active devices in this subcluster is 
ffia (for consistency). Now assuming the serving device to be 
at xq + oq € A^°, the intra-cluster interference at D2D-Rx of 
interest at the origin can be expressed as: 

Tintra = ^ || Xq -f a|l ““ -f ^ Pj || XO + 61 !”“ • 

aeA^O\ao b&B^O 

(32) 

Similarly, the inter-cluster interference can be expressed as: 

Tnter = ^ ^ ^ ||x + a|| “ “t" ^ ^ ||x -f fo|j “j . 

xe(t>c\xo '’S®" 

(33) 

B. Coverage Probability 

We hrst characterize the Laplace transform of inter-cluster 
and intra-cluster interference distributions in the following 
Lemmas. 

Lemma 7. Assuming the serving device to be chosen uni¬ 
formly at random from A^°, the Laplace transform of the 
distribution of intra-cluster interference in (32), conditioned 
on j/Q = ||a:o||, A given by = 

( _ f°° sw~°‘ 

exp - (ma - 1) / fwAwa\vo)dwa. 

V Jo 1 + StUa 

- (Pib) y ^ fwb(.WhWo)dwh^ , (34) 

where fwkiwkWo) = Ricepdf (wfc, j/q; cr^), k e {a,b}. 

Proof See Appendix F. ■ 

Lemma 8. Laplace transform of the distribution of inter¬ 
cluster interference at D2D-Rx of interest in (33) is given by 

'C/inte. (S) = 

/ r°° / / _ su~°‘ 

exp -27rAc / 1-exp -rua / ::/(7a(“ak)dUa 

^ Jo ^ ^ Jo 1 + Slia 

su~°‘ \ w 

-rho ——), (35) 

Jo 1 + SWb ' '' 

where fukiAk\v) = Ricepdf (ufe, z/q; ctI), k S {a,b}. 

Proof See Appendix G. ■ 
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Coverage probability when the content of interest is avail¬ 
able at a device chosen uniformly at random from is given 
by the following Theorem. 

Theorem 4. Using the expression for Laplace transform of the 
intra-cluster interference distribution in (34) and the inter¬ 
cluster interference in (35), the coverage probability at a 
randomly chosen device from the double-variance model is 

roo r(Xi 

P= = / / Ci,^^^^{Pr°‘)Ci,^^^JPr°‘\uo)fR{r\i'o)fvo{i'o)drduo, 

Jo Jo 

(36) 


with, fii{r\uo) =Ricepdi{r,uo-,al) (37) 

where /vo(i^o) = ^ exp(-^). Here, Uk = if D2D-Rx of 
interest is located at denser cluster and = <Jh otherwise. 


Proof The proof follows on the same lines as Theorem 1 with 
slight difference in the distance distributions. By definition, 
coverage probability is 


Pc = EygEfl; 


E 


exp ( Pt (/inter “b /intra)) 


R,Va 


(38) 


where the PDF of serving distance r = ||a:o+ao||, conditioned 
on vq — llxoll, follows Rician distribution [46]. Here, D2D- 
Rx of interest is a randomly chosen device, which is sampled 
from a Gaussian distribution in with scattering variance 
CTa (cb) if D2D-Rx of interest belongs to the denser (sparser) 
cluster. Hence z/g = lls^oll follows Rayleigh distribution. ■ 


It is worth highlighting that the overall performance of 
the double-variance process will depend upon the following 
features: i) serving link distance: it decreases when D2D-Tx 
of interest or D2D-Rx of interest are located in the denser 
subcluster, ii) inter-cluster interference: it decreases with the 
increase of the number of simultaneously active D2D-Txs 
in the denser subcluster compared to the sparser subcluster 
(keeping the total same), and iii) intra-cluster interference: it 
increases with the increase in the number of simultaneously 
active D2D-Txs in the denser subcluster. Here, the first two 
features, i.e., decreasing serving link distance and inter-cluster 
interference, improves coverage probability while increasing 
intra-cluster interference degrades the coverage. 


V. Results and Discussion 
A. Numerical Results 

1) Validation of results: In this section, we validate the 
accuracy of the analytical results, and tightness of the approx¬ 
imations by means of simulations. In all the simulations, the 
locations of cluster centers are a realization of a PPP and 
devices are normally scattered around them. For this setup, 
we set the SIR threshold, /3, as 0 dB, path-loss exponent, a 
as 4, and study the coverage probability for the three cases. 
While the easy-to-compute exact results for the /-Rx and 
baseline cases, given by Theorems 2 and 3, are shown to 
match the simulations exactly, thus validating the analysis, the 
approximations for fc-Tx case given by Corollaries 2 and 3 are 
also shown to be fairly tight. Although the exact expression 
for fc-Tx case, given by Theorem 1, is not straightforward 
to compute numerically, the tightness of approximation given 



Number of the simultaneously active D2D-Txs 


Figure 4: Coverage probability versus number of simultaneously 
active D2D-Txs when a a. — 30, Ac = 50 clusters / km^, Nt = 


Nr = 40. 



Number of the simultaneously active D2D-Txs 

Figure 5: fc-Tx case: Coverage probability versus number of simul¬ 
taneously active D2D-Txs when Ca = 30, Ac = 50 clusters / km^ 
and Nt = 30. 


by Corollary 2 means that it can be used as the proxy for the 
exact result. The results also show that the fc-Tx and f-Rx cases 
lead to higher coverage probability than the baseline case. The 
difference in performance will be even more prominent in the 
ASE that will be discussed in the next subsection. 

2 ) Performance comparison across three cases: Recall that 
there is a clear trade-off between the optimal number of 
simultaneously active D2D-Txs and the resulting interference 
power. While increasing the number of simultaneously active 
transmitters potentially increases ASE, it comes at a price of 
an increased interference. As shown in Fig. 6 for the k-Tx 
case, the optimal number of simultaneously active D2D-Txs 
increases with the decrease in distance from the D2D-Tx of 
interest to the cluster center (i.e., decreasing k). Fig. 5 and 
Fig. 6 show that the coverage probability and ASE are optimum 
when the content of interest for the D2D-Rx of interest is 
available at the closest device to the cluster center. The results 
also show that biasing the content of interest for the D2D-Rx 
towards the cluster center leads to a significant performance 
improvement in both coverage probability and ASE compared 
to the baseline case. On the contrary, it can be seen that both 
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Number of the simultaneously active D2D-Txs 

Figure 6: fc-Tx case: ASE versus number of simultaneously active 
D2D-Txs when = 30, Ac = 50 clusters / km^ and N% = 30. 



Number of the simultaneously active D2D-Txs 

Figure 7: i-Kx case: Coverage probability versus number of simul¬ 
taneously active D2D-Txs when Ca = 30, Ac = 50 clusters / km^ 
and Nr = 30. 

coverage and ASE in fc-Tx case may become worse than the 
baseline case when the content is cached far from the cluster 
center (e.g., fc = 20 case in Fig. 5 and Fig. 6 ). 

Similar trends are observed in Fig. 7 and Fig. 8 for the f-Rx 
case. In particular, the results show that the coverage and ASE 
increase as the distance from D2D-Rx of interest to the cluster 
center is reduced. The impact of £ on the results, especially 
on the coverage probability, is however not as prominent as 
it was in the fc-Tx case. This is due to the fact that while the 
serving distance reduces with decreasing fc, the distances to 
the intra-cluster interfering devices also decrease in general, 
thus leading to a higher intra-cluster interference. 

B. Applications of the Results to Total Hit Probability 

In this section, we use the coverage probability results 
derived in this paper to study the D2D network performance 
in terms of the the total hit probability. We assume that the 
library of popular content for the representative cluster is 
known a priori. It is denoted by the set {ci, C 2 , ...,cj}, where 
the content is ordered in terms of decreasing popularity, which 



Number of the simultaneously active D2D-Txs 

Figure 8: ^-Rx case: ASE versus number of simultaneously active 
D2D-Txs when = 30, Ac = 50 clusters / km^ and Nr = 30. 

means Ci denotes the most popular content. As is usually the 
case, we assume that the content popularity follows Zipf’s 
distribution, i.e., the probability that the content Cj is requested 
by the D2D-Rx of interest is where 7 is 

Zipf’s parameter and J is total number of files [7]. Note that 
the arguments presented in this section are not specific to Zipf 
distribution and can be easily extended to any given popularity 
distribution. The total hit probability can now be defined as the 
probability that the D2D-Rx of interest is able to successfully 
download its content of interest, which in turn depends upon 
two events: (i) this content of interest is available within the 
cluster, and (ii) the D2D-Rx of interest is in the coverage of 
the device that has this content (i.e., SIR > jS). Mathematical 
definition of total hit probability will be provided shortly. 

Due to the limited storage capacities, each device in general 
cannot cache the whole library. For notational simplicity, we 
assume that each device caches exactly one content and the 
number of popular contents is greater than the total number 
of devices, i.e., J > N^. Given that there are transmitting 
devices in each cluster, each device caches one of the N^, most 
popular contents [18]. To evaluate the total hit probability, 
it is possible to have various cache placement and cache 
gathering policies. Due to space limitation, we confine our 
analysis to the following two strategies that directly build 
on the coverage probability results derived in this paper. For 
both these strategies, the D2D-Rx of interest is assumed to be 
chosen uniformly at random from the representative cluster, 
i.e., we confine to the fc-Tx and baseline cases. 

1) Uniform content placement: In this setup, we assume 
that the popular contents are uniformly distributed inside the 
cluster. Recall that the coverage probability when the file of 
interest is available inside the cluster uniformly at random is 
denoted by P® (baseline case). The total hit probability is 

Nt 

Phit=^Pi^,Pc (39) 

i=i 

where P® is the coverage probability given by (28). 































12 



Number of the simultaneously active D2D-Txs 


Figure 9: Total hit probability versus number of simultaneously active 
D2D-Txs when = 30, Ac = 50 clusters / km^, Nt = 30, and 
J = 40. 

2) Cluster-centric content placement: Based on the in¬ 
tuition provided by Lemma 1 for the cluster-centric content 
placement, the most popular content should be placed at the 
transmitting device closest to the cluster center. This means 
that in our setup, ci should be placed at the device closest to 
the cluster center, and C 2 should be placed at the second closest 
device to the cluster center and so on. Hence the caching 
probability of the content cj at the closest transmitting 
device to the cluster center is (j,k < Nt): 

1 j = fc, 

0, otherwise. 

Recall that the coverage probability when the randomly chosen 
D2D-Rx of interest connects to the closest transmitting 
device to the cluster center (fc-Tx case) was denoted by 
Hence, the total hit probability can be expressed as 

Nt 

(40) 

i=i 

where is the coverage probability given by (23). 

We now plot the total hit probability results for the two cases 
in Fig. 9. As expected, the total hit probability is signihcantly 
higher in the cluster-centric content placement case. While the 
shape parameter of Zipf distribution, 7, does not impact the 
results in the uniform content placement case, increasing its 
value improves the hit probability in the cluster-centric content 
placement case. This is because with increasing 7, the most 
popular content is requested more often and since it is cached 
closer to the cluster center, the D2D-Rx of interest connects 
with the devices located closer to the cluster center more often. 
Since the coverage probability in the fc-Tx case for lower 
values of k is significantly higher than for higher values of 
k, this improves the overall hit probability. 

C. Performance of the New Generative Cluster-Centric Model 

In this subsection, we study the coverage probability in 
the double-variance model of Section IV as a function of 
the number of simultaneously active transmitters Wa + rh\j. 



Number of the simultaneously active D2D-Txs (m^ -|- mi,) 

Figure 10: Coverage probability versus number of simultaneously 
active D2D-Txs when and Ac = 50 clusters / km^, CTa = 10, and 
CTb = 30. 

The results are presented in Fig. 10. For each plot, we hx 
either fha or mu and vary the other such that the sum is 
equal to the value on the x-axis. Recall that the analysis for 
this model was performed under the assumption that D2D-Tx 
of interest is sampled uniformly at random from the denser 
subcluster A ^°. It was stated that this assumption will lead to 
a better performance. This can be validated by noticing that 
the coverage probability for the case rfia = 0 (the one where 
both the D2D-Tx and D2D-Rx of interest are in the sparser 
subcluster) is lower than all the cases in which the D2D-Tx 
of interest lies in A^° .Besides, the coverage probability when 
D2D-Rx of interest is located in the sparser subcluster is higher 
than the special case of ma = 0 discussed above, and lower 
than the other extreme in which both the D2D-Tx and D2D- 
Rx of interest are in the denser subcluster (mu = 0). For 
all the cases, we can observe that the coverage probability is 
higher when the number of simultaneously active D2D-Txs in 
the sparser subcluster is higher (keeping Tfia + Wb the same). 
This is because the active D2D-Txs in sparser subcluster cause 
less intra-cluster interference compared to when they lie in the 
denser subcluster. 

VI. Conclusion 

In this paper, we developed a realistic framework for the 
modeling and analysis of cache-enabled D2D networks. By 
modeling the D2D network as a Poisson cluster process, we 
focused on the performance analysis of the cluster-centric 
content placement policies where the content is placed in a 
cluster such that the collective performance of all the devices 
is improved. In particular, we defined and explored following 
two general cases where the location of the D2D-Rx of 
interest or the device that has cached its content of interest 
is parameterized in terms of its location relative to the cluster 
center: (i) k-Tx case: the serving device is the the closest 
device to the cluster center, and (ii) £-Rx case: the receiving 
device is the closest device to the cluster center. Using 
tools from stochastic geometry, we derived the coverage and 
ASE for these cases and compared them with the baseline case 
where both the D2D-Rx of interest and its serving device are 
chosen independently and uniformly at random from the same 





















13 


cluster. The results concretely demonstrated that the network 
performance can be significantly enhanced if either the D2D- 
Tx of interest or the D2D-Rx of interest lie close to the 
cluster center. Based on this observation, we proposed and 
analyzed a new generative double-variance Thomas cluster 
process model that allowed us to tractably capture the fact 
that more intra-cluster interfering devices may lie closer to 
the cluster center. Several system design guidelines for content 
placement, including insights into the effect of realistic content 
placement policies on hit probabilities, have been provided. 

This work has many extensions. From the caching perspec¬ 
tive, it is important to incorporate content popularity distribu¬ 
tion, social interaction between devices and cache constraints, 
such as the memory of the caching devices. From the D2D 
network perspective, it is important to extend the analysis to 
the in-band case where the cellular and D2D transmissions 
share the same spectrum [52]. From the cluster point process 
perspective, it is important to extend the analysis to more 
general classes of cluster processes. From the modeling per¬ 
spective, extensive measurement campaigns must be carried 
out to understand the statistics of real-world clusters, such as 
the ones formed in the coffee shops, libraries, and restaurants. 


Appendix 

A. Proof of Lemma 1 

Recall that the set of possible transmitting devices inside 
a representative cluster is denoted by G {l,2,...,iVt}. 
Assuming the file of interest is located at xq -f Sfc, we have 


k* =arg maxj,g_^-o E[l{SIR(||a:o + Sk\\) > P}] 


(a) „ 

= arg maxj,g_;^xo E 

= arg max^gjy-^^o E 


hao\\xo-\-Sk\\ 

.It - KoWxq + Sfe||"“ 

iko+5fcir° 

j^It - lla^o + Sfc||““ 


= arg maxj,g_^-o E [Xk > /3] 


where I-r in (a) is the total received power at the D2D-Rx 
of interest from all the transmitters in the network, and Xk 
is defined as the received SIR from the closest serving 
device for the ease of notation. Since It is not the function of 
k, Xk >st Xj whenever ||a;o + Sj|| >st lla^o + s/c||, where >st 
denotes first order stochastic dominance (or usual stochastic 
order). Note that since xq is sampled from zero mean complex 
Gaussian random variable in the density function of 
r = IIxq + Sfell conditioned on tk = ||sfe|| follows Rician 
distribution with CDF Fji{r\tk) = 1 — Qi{^, f-), where 
Qi{a,P) is the Marcum Q-function defined as Qi(a,P) = 
ye~^^ Io{o(y)dy. It turns out that Qi{a,P) is mono- 
tonically increasing in a [53, Property 11], which implies 
Qi(^, is monotonically increasing in tk, which implies 
R{tk) <st R{tk + 1). This implies ||a;o + Sj|| >st H^o + si|| 
Vj f 1, which completes the proof. 


B. Proof of Lemma 2 

Let the angle between the intra-cluster interfering device 
and the D2D-Rx of interest, as seen from the cluster center. 


be 6. While the angle is uniformly distributed between [0, 27r], 
the “direction” is not important for the distance calculation, 
which means it suffices to assume it uniformly distributed 
between [ 0 , 7 r]. Now, the CDF of the distance between these 
two devices, conditioned on vq = ||a:o|i and t = ||a|| is 


Fwiw\iyo, t) = F[W < w\uo, t] (41) 

P [j/q -P — 2iyot cos 9 < w‘^\vo,t\ 
cos9 > ^^ 


2vQt 


ib) 


< cos 


-1 


^ f^ — 


F) 1 -1 

= — cos 

TT 


2vot 

2vnt 




\vo,t 


where (a) follows from the cosine law, (b) follows from the 
fact that cos“^ is monotonically decreasing function, and (c) 
follows from the fact that 9 ~ Unif[0, tt]. Now, the conditional 
PDF of fw{w\vQ,t) is obtained by differentiating over w as 
follows: 


fw{w\uo,f) 


1 

TT 


wjv^t 



ko 


t\ < W < Vq -\-1. 


where |.| denotes absolute value. Using the fact that devices 
are normally scattered around cluster-center, the distance from 
intra-cluster devices to the cluster center, i.e., t = ||a|| is 
Rayleigh distributed with parameter CTa, which implies that 
the PDF of distance fi„ of an intra-cluster interferer in the set 
i.e., fin < tk is truncated Rayleigh distribution 


/Ti„(Ln|ifc) = /t (finikin < 4) 


friUn) 
Fxitk) ’ 


< tk, 


where /t( ) and FV(.) are the PDF and CDF of Rayleigh 
distribution with parameter da respectively. Similarly, the PDF 
of distance fout, where tout > ffc, is 


f Tjut {t out \tk) = frit out \Tout > tk) 


/T(fout) 

1 — Fritk) 


^out 


> tk- 


Using the same approach as [46, Lemma 4], the conditional 
i.i.d property of w = ||a;o + a||, conditioned on vq = ||a;o|| and 
t = ||a|| can be formally proved. 


C. Proof of Lemma 3 

Assuming that the file of interest is located at the k*^ closest 
transmitting device to the cluster center, we divide the set of 
simultaneously active intra-cluster devices into three subsets: 

= {^fj[’,ao,-4nnt}. Here oq denotes relative location 
of the serving device to the cluster center with distance 
||ao|| = tk away from it, where Aff = {a|||a|j < tk}, and 
-^out = {®| ll®ll > tk} are the set of devices closer and further 
than serving device from cluster center, respectively. For this 
setup, the Laplace transform of distribution of intra-cluster 
interference 

='E|^exp(^-s ^ ||a:o + a|r“)] 

aGA.^0 \aQ 
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= E|exp(^-s(^ ha^Jxo + a\\~ 

+ Y1 ^a,oll®o + “ll~“))] 

“e-^oSt 

= e[ exp(-s/ia,J|a:o + a|r“) 

aeAf^ 

X n ®xp(-s/ia,J|a;o+ a|r“)] 


(42) 


(43) 


= E 


n 


^ n 

1 + s||a;o + a||““ ■*'L ^ 

aeA^° aeA^° 


1 

+ s||a;o + a|| 


is exponential distributed with mean unity, (c) follows 
from expectation over number of devices that is Poisson 
distributed conditioned on the total being less than TVt — 1 
along with the fact that locations of devices conditioned on the 
location of cluster center, xq, are i.i.d, and (d) follows from 
converting Cartesian to polar coordinates using the fact that 
fw the density function of distance from interfering 
devices to the D2D-Rx of interest conditioned on i^q = ||a;o||, 
has Rician distribution given by (9). Now under the assumption 
fn-a ‘C iVt, the Laplace transform of intra-cluster interference 
distribution can be approximated as: 


(c) 


iVt-l g„ 


n=0 1=0 


Wo 1 + 


—/w ( w I t'o, Ln)/Ti„ (Ln I i fe )dwdii, 




‘ ~ ^ exp - (Wa - 1) J 


— fw{w\i^o)dw). 


1 -I- sw~°‘ 


Mi^(s,tkWo) 

fw (ro| r'o, fout (^Tout (tout I (dwdiout 


(7^(1-pr 


7(1 -p;n- ffm, 1 -fflrn)^ 

¥{L = l\L<gm) 


-^out )^fc I ^o) 

(ma - 


P(iV = n|iV<iVt-l) 


with w\t = \v — t\. 


w}i = r-o + (ii 


E. Proof of Corollary 3 

n — l 

Under Assumption 1, the correlation corresponding to the 
common distance t'o = |lxo|| is ignored and hence the Laplace 
transform of the interference distribution can be approximated 
by Corollary 1. Furthermore, the density function of serving 
distance r = |jxo -f Sk|| € ffi+, needs to be evaluated only 
conditioned on the distance of D2D-Tx of interest (i.e., 
closest device) to the cluster center tk = ||sk||- Now using the 
fact that Xq e ffi+ is zero mean complex Gaussian random 
= |jyp _ variable, the density function of z = xq -f Sk € where 


= 1^0 + tout, P = tEa’ 5m = min(n,fc - 1), and z = (^ 1 ,^ 2 ) conditioned on Sk = (sfci,SfcJ (where tk = 


. Here (a) follows from definition 


■^out ~ ‘'U I ‘-outj F — 

0 _ —1 (ma— 

S — 2^j=0 JT 

of Laplace transform, (b) from the fact that ~ exp(l), 
and (c) from converting Cartesian to polar coordinates by using 
distance distribution given by Lemma 2 along with conditional 
i.i.d. property of fwiw\iyo,t). Then, the result follows by 
expectation over number of devices, where the number of 
devices closer than serving device to the cluster center, i.e., 
I, is binomial distributed conditioned on the total being less 
than (/m = min(n, fc — 1). This condition is due to the fact 
that I is always smaller than k — 1 (since serving device is 
located at device in fc-Tx case) and total number of active 
devices, i.e., n, where n is Poisson distributed conditioned on 
total being less than Nt — 1. 

D. Proof of Lemma 4 

The Laplace transform of distribution of intra-cluster inter¬ 
ference can be derived as: = E [exp (—s/intra)] 


) can be expressed as: 


’fc2 


/z(zi,Z2|Sk) = 


(Zl-Sfcj2 {Z2-Sk^f 

2cr2 2cr2 


Since, we are interested on distribution of r = |jz||, we define 
zi = rsin5, and Z 2 = rcosO, where 9 = arctan(|^). Now, 
Jacobian matrix is used to convert Cartesian coordinates to 
polar coordinates as follows: 


/it,e(L^'|sk) =/z(2i,z2|sk) X 


d 


Zl,Z2 


r. 


(44) 


where d 


/ \ 

dzi 

dzi 

/ 21,22 \ _ 

7 J “ 

dr 

OZ 2 


dr 



= r, and hence joint 


distribution of (7?, 0) is /K,e(r’, 5|sk) = 

(r cos 9- SkJ"^ (r sin 9 - Sk^)‘^ 


= E 


[exp 
=E^i^o [ 


-s E 

a^A^O \ao 


^h„ 


aGA^O \aQ 


_J|a:o + a|| 

[exp {-sha^, 


0] 


27rcr2 


exp - 


= ^exp - 


Ifo -I- a|| 




-I- t 


2(t2 


2^2 

1 

27r 


exp 


2o-2 

rski cos 9 + rsk 2 sin 9 


Therefore, the conditional marginal distribution of R is 


‘=Wo[ n 


1 


aGA^O \ao 


1 -I- s||a;o + a|| 




fR{r\tk) = / fR,e{r,9\si^)d9 = — exp - 


2a2 


(^E 

n=0 

Nt-1 

(JE 

n=0 


]g2 1 -h s||xo + Cl 
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-fA{a)da^ 


fw iw\uo)dwJ 


withe = E":; 


1 + sw~°‘ 

Nt-1 


■ (ma - 

"(rna - 


c2-k 




/ rskt ^ 9 

I 


d9, 


j=0 


, where (a) follows from the 


where conditioning on tk = ||sk|l instead of Sk suffices. The 
rest of the proof follows on the same line as the proof of 


definition of Laplace transform, (b) follows from the fact that Theorem 1. 



15 


F. Proof of Lemma 7 

Since the two sets A^° and are independent, the 
Laplace transform of intra-cluster interference distribution, 
'^hntra('Sl^'o) = E [cxp (—sJintra)], Can be derived as follows: 


= E|^exp(^-s ^ ha^o\\xo + a\\ 

aGA.^0 \aQ 

+ ^ hb^g\\xQ + b\\ 

beB^o 

=E^-o [ Eh„^^ [exp {-sha^g 11*0 -f a|r“) ]] 


a£A^0 \aQ 


(c). 


X [ n [expHVJ|*o+6|| “)]] 
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E^ 0 [ P _|_fl,||-a] 
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■ Eb=:o 
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+ s||xo -f ( 

1 
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(45) 
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exp - (ma - 1) 


o -I- *0 


- mb 


/k2 


r 2 1 -f s||a;o -I- a|| 
+ “ /sWdfe) 


-fA{a)daj 
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exp ( - (ma - 1) / 7 /Wa(Wa|l'o)d«la 

V Jo 1 + SWa 


- mb 


1 -I- SWh 


1 -I- SWa 
- fwkiwbWoHwbJ 


where (a) follows the definition of Laplace transform, (b) 
follows from the fact that the interference from set and 
are independent, (c) follows from the fact that 
and are exponential random variables with mean unity, 
(d) follows from probability generating functional (PGFL) 
of Poisson distribution where /b( 6) = exp ^ ~J |^2 ^ 

and (e) follows from converting from Cartesian to polar 
coordinates and some algebraic manipulation similar to the 
derivation of Rician distribution in the proof of Corollary 3. 


G. Proof of Lemma 8 

Laplace transform of the distribution intra-cluster interfer- 


ence at D2D 

-Rx of 

interest Ci- 

ntor ('®) 

is 


=e[ exp (—s 

E 

\ E 

\x + a|| 


E ^bAx + b\\ 



re</)c\a;o 


bGB^ 

= 1 

n 

Eyix 

[ n 

[exp(- 
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L n 1 -(- s 
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= exp ( - 2-7rAc / ( 1 - exp ( - rUa / --(Ma|!^)dUa 

5 Jo ^ ^ Jo 1 + S«a 


- mb 


1 -I- SMb 


1 -I- SUa. 

Aubi'^b\v)dub'juduj'j 


where (a) follows from the fact that A^ and B^ conditioned 
on the location of cluster centers {a;} are independent, (b) 
follows from the fact that ha^ and hb^ are independent 
exponential random variables with mean unity, and (c) follows 
from the PGFL of Poisson distribution. Note that Lemma 5 is 
a special case of Lemma 8 when fhh = 0. 
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